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Abstract 

We investigate a model of dynamical electroweak symmetry breaking via a dual 
gravitational description. The gravity dual is obtained by embedding a D7-D7 pair 
of branes into a type IIB background that is dual to a walking gauge theory. We 
develop further a previous study of this model. In particular, we show that there is a 
nontrivial relation that needs to be satisfied in order for axial-vector modes to exist. 
Furthermore, we compute explicitly the electroweak S parameter. The result is 
positive-definite and, as was to be expected, much smaller than in earlier QCD-like 
D-brane constructions. We also find the masses and decay constants of the vector 
and axial-vector mesons in this model. This allows us to obtain another estimate 
for S by summing the contributions of the discrete states. It is noteworthy that, 
in contrast to previous holographic studies, the sum of the first several lowest-lying 
states does give a very good approximation to the full answer. 
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1 Introduction 

The spontaneous breaking of electroweak symmetry and the ensuing generation of gauge 
boson masses are key ingredients of modern particle physics. Uncovering the origin of 
this electroweak symmetry breaking (EWSB) is one of the most important goals of fu- 
ture experiments. In the Standard Model the mass generation is due to the presence of 
an SU{2) doublet Higgs field, comprised of four real components. Spontaneous EWSB 
occurs when this field acquires a non-vanishing vacuum expectation value. Three of its 
components then give rise to the would-be Nambu-Goldstone bosons of the spontaneously 
broken symmetry, which become the longitudinal third polarizations of the massive 
and Z bosons. The remaining scalar degree of freedom is the Higgs boson. 

Notably, however, the Higgs boson has not yet been discovered experimentally. Fur- 
thermore, as it is a fundamental scalar, its mass diverges due to quantum corrections. 
This problem can be significantly alleviated in supersymmetric extensions of the Stan- 
dard Model. Namely, supersymmetry (SUSY) tames the quantum corrections to the Higgs 
mass by providing automatic cancellations between bosonic and fermionic superpartner 
contributions. Despite that however, an initial fine-tuning is still necessary. This lack of 
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understanding of the hierarchy between the electroweak (Higgs mass) scale and the cutoff 
(Planck) scale is known as unnaturalness problem. In addition, the question of SUSY 
breaking remains to be settled. 

All those reasons have motivated the search for a more conceptually appealing alter- 
native. A candidate for such an alternative, that leads to EWSB without any fine-tuning, 
is offered by technicolor (TC) theories [1]. In the latter, instead of the Higgs mechanism, 
one uses dynamical chiral symmetry breaking in a new strongly-coupled gauge sector. 
Initially, TC models were thought to have dynamics similar to QCD, only at a higher 
energy scale. However, such TC theories face important problems when compared to phe- 
nomenology. Namely, they predict too large a value for the Peskin-Takeuchi S-parameter 
[2]. Furthermore, the introduction of extended technicolor (ETC) interactions [Hj, in order 
to generate fermion masses, leads to unacceptably large flavor- changing neutral current 
amplitudes. 

A phenomenologically promising modification of the above ideas is provided by walking 
TC (WTC) theories A characteristic feature of the latter is that the gauge coupling 
runs much slower than that in QCD between the TC and ETC scales. This is due to 
a different field content compared to that in QCD-like models. The presence of such 
a walking (almost conformal) intermediate region is expected to help solve the problem 
with the magnitude of the flavor-changing neutral currents. Starting from [3], vector 
gauge theories, that exhibit strong IR fixed point behavior, have been used to build 
WTC models]^ It has also been conjectured that WTC theories should predict smaller 
S-parameter than QCD-like TC ones, due to the markedly different bound state spectra 
O [8]. Recall that, in QCD-like technicolor the S-parameter was obtained by scaling-up 
to the TC scale the measured spectral functions of QCD. In WTC, however, such an 
estimate is clearly incorrect. Instead, one has to address the computation anew, while 
working at strong coupling. This is, obviously, beyond the realm of perturbative QFT 
methods and thus is the main issue that has long hindered the investigation of walking 
technicolor. For reviews of WTC, see 

Modern developments in gauge/gravity duality provide exactly the tools to address 
the above problem, namely to study gauge theories at strong coupling. What we mean 
by this is the following. Gauge theories, that can be realized as the worldvolume degrees 

"'in 6_ the techniquarks were placed in higher representations, in order to achieve walking with a 
smaller technifiavor symmetry group and thus to reduce the value of S. 

^We should note that in recent years nonperturbative methods, like lattice simulations [TD] and solving 
bound state equations |TT], have been utilized to estimate the spectra of non-QCD-like theories. These 
investigations show reduction in S compared to QCD-like models. 
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of freedom of D-brane configurations in string theory, have dual descriptions in terms of 
the gravitational backgrounds created by those branes. Most importantly, the strongly- 
coupled regime of the field theory corresponds to the weakly-coupled regime of the gravity 
dual. Hence, one can use classical supergravity to obtain answers to questions about 
strongly-coupled gauge theory. This remarkable duality has been utilized to investigate a 
wide-ranging set of problems, from hydrodynamics [12] to superconductivity [I3]|^ Most 
relevant for us, [131 US] built a gravitational dual of a gauge theory that exhibits many 
of the characteristic features of QCD. In this holographic QCD model, chiral symmetry 
is spontaneously broken due to a particular U-shaped embedding of flavour D8-D8 pairs 
into a background sourced by color D4-branes. 

The above Sakai-Sugimoto model can easily be adopted to give a gravitational dual of 
QCD-like technicolor. The latter, as well as other similar models obtained from various 
brane configurations in type IIA/B string theory, has been studied in [16] |^ With this 
method one can compute explicitly the Peskin-Takeuchi S, T and U parameters |2J , that 
are important electroweak observables on a par with the fine structure constant and the 
mass of the Z-boson. As is to be expected though, the QCD-like technicolor models of 
[16] are incompatible with present observational bounds. In [18], a gravitational dual of a 
walking technicolor model was considered. A key ingredient was the type IIB background 
of [in], which is dual to a walking gauge theory. One can show that there is a U-shaped 
embedding, a la Sakai-Sugimoto, of a probe D7-D7 pair in that background [IB]. Thus, 
one can realize geometrically chiral symmetry breaking, which then induces electroweak 
symmetry breaking. 

Here we will study further this gravity dual of walking techinicolor. First, we will show 
that there is a non-trivial constraint between the three parameters of the model, which has 
to be satisfied in order for regular axial-vector modes to exist. Hence, the number of free 
parameters is with one less than initially expected. Then we will compute analytically 
the S parameter via the exact expression in terms of non-normalizable modes derived 
in [ini US]- Unlike [IS] however, our present method for finding the non-normalizable 
solutions will not leave any undetermined constants. The explicit answer we obtain for 
S is positive-definite and compatible with the current experimental bounds, although it 
does exhibit an unexpected functional dependence on the length of the walking region. 

®It should be noted though that there is a large variety of AdS /CFT- inspired phenonienological models, 
whose string theory embedding is an open problem. 

^Note that there is a huge, (loosely) inspired by AdS/CFT, literature on holographic technicolor 
models |17j , which do not have a consistent string theory embedding. 
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We also compute the masses and decay constants for the spectrum of vector and axial- 
vector mesons in this model. This allows us to obtain yet another estimate for the value 
of the S parameter, by summing the contributions of the discrete states. The result 
agrees completely with the answer from the non-normalizable modes, in contrast to the 
QCD-like technicolor models of [16]. In fact, the sum over discrete states converges rather 
rapidly. Hence, summing several lowest lying states does give a very good approximation 
to the full answer. Finally, it is worth noting that we can relate directly to each other, 
i.e. without any of the model parameters present, the physical observables S, F^^ and mj^, 
where F^^ is the technipion decay constant and is the mass of the n^^ vector state. 
This relation implies, in particular, that for a single techniflavour pair the lightest vector 
meson mass m\ > 5 TeV, in order for the experimental bound on the S parameter to be 
satisfied. 

2 Gravitational dual 

One can build gravitational duals of technicolor models in the following manner. Let us 
consider a background sourced by a stack of Ntc technicolor branes in type IIA/B string 
theory. To add techniflavour degrees of freedom, we need to include an additional set 
of NxF probe branes and N^p probe antibranes. If there is a U-shape embedding a la 
Sakai-Sugimoto ^1] of the techniflavour probes into the technicolor backrgound, such that 
the probe branes and antibranes merge at some finite radial distance, then we have chiral 
symmetry breaking. This is because the techniflavour group Ul{Ntf) x UrINtf), arising 
from the separate branes and antibranes, gets broken to the diagonal subgroup U{Ntf) 
due to the merging. This chiral symmetry breaking then induces electroweak symmetry 
breaking via an appropriate choice of embedding of the electroweak SU{2) x f/(l) group 
into the techniflavour Ul{Ntf) x Ur{Ntf)- The above general set up has been utihzed to 
study a variety of QCD-like technicolor models [16]. In [18] it was applied to build a model 
of walking technicolor. An essential ingredient in that regard is the recently discovered 
gravitational background dual to walking behaviour [IS]|^ In this section we review briefly 
the background of [I9] and then discuss in more detail the U-shaped techniflavour probe 
embedding studied in [18]. Finally, we show that there is a nontrivial condition, that has 
to be satisfied in order for regular axial-vector modes to exist in this model. 

^We should note that this is within the class of the solutions studied in [20] , 
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2.1 Walking technicolor background 

The background of [191 is an A/" = 1 solution of the type IIB equations of motion with non- 
trivial RR 3-form flux and constant dilaton. It is a deformation of the Maldacena-Nunez 
background |2T], that still arises from D5-branes wrapping an S*^. The ten-dimensional 
metric is: 



A 



dxi 3 + 



+ 



cPi{p) coth(2p) 
4 



wher43 

A = 



sin^ a 



1/4 



d^l + d^ll + 



cosh(2p) 



Piip) 



dp 



Pi{p) = (cos'^ a + sin^ a (sinh(4p) 



(2.1) 



(2.2) 



cui = cos ijjdO + sin ijj sin Odcp , 
u}2 = — sin ipd6 + cos ip sin 
0)3 = dip + cos 9dip , 



ui = d6 , 
UJ2 = sin 6d(p 
Us = cos 9d(p 



(2.3) 



and 



ddl = cof + col ? 



d9^ + sin^ 9dip^ 



(2.4) 



Here the constants c and a are parameters characterising the solution and, as in [19], we 
have set a' = 1 = Qs- Also, we will not write down the explicit form of the RR 3-form 
flux, since it will not be needed for our purposes. 

The walking region corresponds to the interval p G (l,p*), where the upper end is 
approximately given bj|^ 



p*'^ - (log 2 + 3 log cot a) 



(2.5) 



This can be seen by investigating the dependence of the technicolor gauge coupling 
on the radial variable p, which can be read off from the DBI action of a supersjmametric 
D5-brane probe in the above background |19j : 



arcsinhl . . ,„ . 

g V sinh(2p) y 



cPi 



c (cos^ a + sin^ a (sinh(4p) — 4p)) 



1/3 



(2.6) 



^In the expression for A we have corrected a typo in [19]. 
We will comment more on the value of in Section jsj 
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Let us also note that in the walking region the metric (2.1) simplifies to: 

c cos OL ( tan^ a e^^ 



dxi 3 + 



4 V 



(4dp^ + (wg + wg)^) + dnl + dn 



(2.7) 



Finally, the validity of the approximation, in which the above solution has been found, 
requires that the parameter c ^ 1, whereas for the existence of a walking region one 
needs that (3 = sin^ a <iC 1. 



2.2 Techniflavour probe branes 

To add technifiavour degrees of freedom, we will consider probe D7 branes in the back- 



ground (2.1). In [18] it was shown that one can find a probe embedding of the necessary 



U-shape type for the following ansatz: The D7 worldvolume directions are taken to be 
the four space-time dimensions x^, the radial variable p and the 3-sphere parameterized 
by ip, 9 and (p. Furthermore, the transverse coordinates 9 and ip are functions of p only, 
i.e. 9 = 9{p) and ip = 

The embedding functions 9{p) and (p{p) are determined by the field equations that 
follow from the DBI action for the probe branes or, in other words, from the Lagrangian 



C = const ^ — det(5'8 



(2. 



where is the metric induced on the D7-brane world-volume, i.e. the metric (2.1) with 
9 = 9{p) and ip = ip{p) substituted. More explicitly, we have 



V-detigs) = [m + gip) {9^ + sin^ ^^J)] 



2M V2 



(2.9) 



where 9p = d9/dp , 'Pp = dip /dp and 

1 



m 

9{P) 



256 
1 



X A^c*Ptip)P[\p) coth'(2p) 
''"Vi'(p)P{(p)coth(2p) . 



(2.10) 



256 2 

It is easy to see that the equation of motion for 9{p) is solved by 6' = tt /2 for every 



p, just like in the simplified walking metric (2.7) that was studied in [18j. Then the 



Euler-Lagrange equation for ip{p), that follows from the Lagrangian (2.9), is: 

d dC 



dp dipp 



0. 



(2.11) 



'^This is very similar to the embedding of D7 probes in the conifold, considered in ^22j. 
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which imphes that: 



dC 

dipp 



C 



with C being a constant. More exphcitly, (2.12) has the form 



just an algebraic equation for (pp that can easily be solved: 



(2.12) 
C, which is 

(2.13) 



' 9{P) [9{P) - C'] ' 

To fix the constant C, let us look at how things simplify in the walking region. Due 



to the simplified metric (2.7), the Lagrangian acquires the form: 

C = const X 62^(456^" + iflY^^ , 

where B 



tan^ a 



Hence, the general equation (2.12) reduces to 



implies that 



,1/2 



(fp = ±2\fBC- 



1 



C2\l/2 • 



(2.14) 
C, which 

(2.15) 



Comparing this with the walking region embedding equation (4.19) of [iSj, we conclude 
that: 



C 



(2.16) 



where po is the value of p at which the D7 and anti-D7 branes merge, i.e. the tip of the 
U-shape embedding. To recapitulate, the embedding function p{p^ is determined by 



3 4po 



/(p) 



(2.17) 



g{p) \g[p) - e^po] " 

So far, we have only addressed the embedding of the techniflavour branes into the 



background (2.1). In order to study vector and axial- vector mesons, that are necessary 



ingredients in a technicolor model, we need to include the fluctuations of the worldvolume 
(techniflavour) field strength Fab of the D7 branes, with a, 6 = 0, 1, 7. In other words, 
we have to consider the DBI action 

Sdi = -T J dp dn^ e- V- det((?,, + 2na'Fab) , (2.18) 

where Qab is the metric induced on the D7 worldvolume, Q3 is the 3-sphere wrapped by 
the D7 branes and is the dilaton. To leading, quadratic, order in Fab the above action 
becomes: 



S 



D7 



-T{27ra'y J d^x dp dn^ e" det(^8) Tr g^'g^^^FacF^ • (2.19) 



Using the embedding solution (2.17), together with 6 = n/2, and integrating over ^3, we 
obtain: 

Sd7 = -lj d'^dp [a{p)F^,F^^- + 2h{p)F,pF%] , (2.20) 

where p,u = 0, 1, 2, 3 and n = "^(^^^ ) ^3 with V3 being the volume of the compact cycle 
wrapped by the D7 probe brane. Furthermore, the functions a{p) and b{p) are: 

Kp) = ^Af^P)+9{p)v'^"'9'' . (2.21) 
where the inverse world-volume metric component (/^^ is given by 

^ A^c^A^P; coth^(2p) ^ 2 

128 (/ + g^f) Ac [P[ + \P, tanh(2p) y,^] " ^ " ^ 



Finally, let us determine explicitly the coefficient k in (2.20). Since V3 is the surface 
"area" of a unit 3-sphere, it is given by V3 = 27r^. Then, recalling that we work in units 
of = a' = 1 = , we have: 



(2vr)7/f ^ 2(2 



TT 



To conclude the description of the gravity dual set-up, let us also recall that the 
vector and axial-vector modes we are interested in arise from the following expansion of 
the techniflavour gauge potential A^(x,p): 

A,{q,p) = V,(g)^° (g^p) + A,{q)rA{q\p) + E Ki<l)^yr^ip) + ^»^aM , (2-24) 

n 

where we have Fourier transformed in x^. Furthermore, the modes in the sum over 
n are the normalizable ones, while tpy and are the non-normalizable modes that 
correspond to sources for the vector and axial-vector field theory currents, respectively!^ 
An important point is that, as in the Sakai-Sugimoto model [H], the vector modes are 
those that are symmetric w.r.t. refiection (on the probe brane embedding) around po, 
whereas the axial modes are those that are antisymmetric w.r.t. this operation. For more 



i^Note that, just like [12], in (2.24 1 we neglect an additional term piopoitional to the derivative of 
the technipion field, as that is not important for our purposes. Such a term arises from the gauge 
transformation that sets Ap = Q. It plays an essential role in obtaining the technipion's kinetic term and 
its interactions with other fields in the four-dimensional Lagrangian, just as in [15]. 
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details on the expansion modes see [18]. Here we will only write down, for future use, the 
normalization condition 

K / _dp a{p)'4)n'4'm = Km , (2.25) 
JD7+D7 

that is needed (both for V and for A modes) in order for the four-dimensional action to 
be canonically normalized, as well as the equations of motion 

dp[b{p)d,Mp)] = -mla{p)Mp) (2.26) 

and 

d,[bip)d,K^{q',p)] = -q^a{p)i^%{q\p) , (2.27) 



that follow from the action (2.20) with (2.24) substituted. 



2.3 Condition for regular axial- vector modes 

Let us now take a closer look at the equation 

h{p) r{E\ p) + h\p) i^\E\p) + E\{p) i;iE',p) = 0, (2.28) 

where ' = dp and E denotes any of mv'„,m^,^,g. We will show that there is an in- 
teresting condition, that has to be satisfied in order for this equation to have regular 
axial- vector solutions ipA„- In order to do that, it will be convenient to change variables 
from p to a worldvolume coordinate z that distinguishes between the D7 and the D7 
branches. Furthermore, since vector (axial-vector) modes are defined via their symmetry 
(anti-symmetry) with respect to reflection around the point p = po on the D7-D7 embed- 
ding, the coordinate z should be such that p = Po corresponds to 2; = 0. A suitable choice, 
discussed in [TS], is = p^ — pg, where z > (z < 0) gives the D7 (D7) branch. Another 
choice, that will turn out to be natural in our context, will be introduced in Section |4j For 
any change of variable p p{z) though, one finds that the equation of interest becomes: 

b{z) ip'^E"^, z) + y^z) i)'{E?, z) + E^a(z) i){E? , -2) = , (2.29) 

where now ' = dz and ^18] : 

a{z) = a{p)dzp , K^) = |^. (2.30) 

Ozp 

Note that, clearly, physical observables should not depend on the choice of z. This is 



ensured by the relations (2.30), as will become evident in Section 3 On conceptual level 
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the independence of physical results on the concrete function p{z) is due to the fact that 
different choices of z differ only by a worldvolume coordinate transformation. 



Now, we are looking for two types of solutions of (2.29): vector modes, which are 
defined by being symmetric under z — )■ —z and thus satisfy ipy{0) = 0, and axial-vector 
modes, that are antisymmetric under z — )■ —z and so ^/'a(O) = 0. As it turns out, there is 
a nontrivial condition for the existence of regular antisymmetric solutions. This is due to 
the different behaviours at 2; = of the functions of interest when g{po) 7^ e^^° and when 
g{po) = To see this, let us consider in more detail the small z expansions of the 

various relevant expressions. We will only need that dzP ~ 2; at small z for any change of 
variables of the form z"^ 
regular at p 



f{p — po). It is easy to realize that /(p) and g{p) in (2.10) are 
Po, or equivalently at 2 ~ 0. On the other hand, in (2.17) behaves as: 

const + C(2^) , when g{po) ^ e^''° 



+ ... 



when g{po) = e^''° 



2.31) 



Using this in (2.21 ) and (2.22), one can extract the small z behaviour of a{z) and b{z). As 
a result, one finds that there are always even (vector) solutions of (2.29), that are regular 
at z = 0. However, regular odd (axial- vector) solutions exist only for g{po) = e^^^ . When 
this condition is not satisfied, the odd modes develop a logaritmic singularity at 2; = 0. 

Given the importance of the condition g{po) = e'^^° for the existence of a technicolor 
model, let us write it down more explicitly: 



256 



[8(l-/32/3)3/2+4/3(e 



4po 



-4po 



5P0 



2c^l3e 



4po 



(2.32) 



where for convenience we have introduced the notation (3 = sin^ a and used that cos^ a = 
(1 — /3^/^)'^/^. This relation allows us to solve for one of the three parameters /3, c and 
Po in terms of the other two. So we end up with one less free parameter than originally 
expected. Let us also comment on the consistency of this condition with the necessary 
parameter ranges, more precisely: 



c » 1 



and 



13 « I, 



(2.33) 



as recalled under (2.7). It is easy to see, keeping po fixed for simplicity, that if one 



increases the value of c, then (2.32) leads to decreasing the value of /3. So there are 



infinitely many solutions consistent with (2.33). This is even more evident upon realizing 



that the constraint (2.32) simplifies significantly at leading order in small (3, due to the 
fact that Po is always at least of order 1. Namely, we have: 

yS^lGc^. (2.34) 
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3 S-parameter: exact result 



In this section we will compute the electroweak S parameter using the exact formula 
of HE], that follows from the gauge/gravity duality statement that the generating 
functional of the field theory correlators is the dual gravitational action. This relation 
allows one to extract a compact answer in terms of the non-normalizable modes introduced 
in the previous section. To find those modes, we will use a different approach compared 
to [18], which will not leave any undetermined constants at leading order in small /3. 



3.1 Preliminaries 



The S-parameter is defined as [2]: 



d 



s = -4.-(n.-n,) 



(3.1) 



where 11^ and 11^ are the two-point correlators for the vector and axial-vector currents 
respectively. In [18J, it was shown that: 



2k [fe(p)^°(g',P)5p^y(?',P)] 

2k [bip)rAiq',p)dMq',p)] 



p=oo 



p=oo 



(3.2) 



where the constant k is given in (2.23), ipyj^ are the solutions of (2.27) and, finally 



b{p) is the same as in (2.21). Recalling that the technipion decay constant is given by 
= 11^1(0), we can see from (3.2) that: 



2k [b{p)r^{0,p)d.^l{0,p)] 



p=oo 



(3.3) 



However, applying the above formulae to the present case requires a slight modification 
(something that was not realized in [18] )• The reason is that the background of interest has 
a walking region and a UV region and the transition between these two regions is related 
to a change of spectrum in the field theory. That was not the case for the duals of regular 
(i.e. QCD-like) technicolor, for which there was just a single region and so the integration 
in p was going all the way to 00. In the present case however, such an integration would 
not encode correctly the relevant physics, since S should be determined by the modes in 
the walking region; the modes in the UV region, on the other hand, should somehow be 
related to extended technicolor So in the following we will introduce a physical cut-off 



^•^Whether (or not) the background of [19], as well as the necessary techniflavour probe embeddings, 
require some modification in order to properly account for this extended field theory is an interesting 
open question, that we hope to address in the future. 
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Pa, which is the upper end of the walking region. Hence, instead of (3.2) and (3.3), we 
should compute the S-parameter and from: 



S = -Stt/s: b{p) ^ [ry{q\p)d,rv{q\p)-^W,P)d,rAiq\p)] 



p=PA,q'^=o 



and 



F^ = 2k [6(p)^°(0,p)a,^°(0,p)] 



P=PA 



(3.4) 



(3.5) 



respectively. Note that this is similar to [23j, where a physical cut-off was introduced in 
the computation of glueball and meson spectra in (deformations of) the Maldacena-Nunez 
background, that is dual to A/" = 1 SYM, 
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In view of the above, in order to satisfy the boundary condition ^^^{pa) = 1 we have 
to take: 

To incorporate this normalization automatically, we write Ti^ y as: 



2Kh{p)dp\og[ijy{q\p)] 
2Kh{p)d,\og[^l,A{q\p)]\ 



P=PA 



P=PA 



(3.7) 



From now on, we will use these expressions to evaluate (3.1). Before turning to the 



computation of the S-parameter though, let us first make a few comments about the 
cutoff Pa- 



In [19] it was shown that the upper end of the walking region is roughly given by (2.5 ). 



This estimate is obtained by investigating (2.6 ), i.e. the behaviour of the technicolor gauge 
coupling g"^ as a function of the radial variable p. One can see that the walking region 
corresponds to the sum in the denominator being dominated by cos^ a. Hence, a rough 
estimate for pA, the upper end of the walking region, follows from the condition that 
sin^ a (sinh(4p) — 4p) becomes comparable with cos^ a. In other words: \e 



cot^ a, 



since p is always at least of order 1 or greater; this is the condition (2.5). However, this 



cutoff is, in fact, in the middle of the smooth transition between the walking and the UV 
regions. A more careful estimate, that removes this transition region, is provided by 

1 



--e^p^ = 0.02 X cot^ a . 



(3.8) 



^''As an interesting aside, let us also mention that the introduction of a cut-ofF in the QCD-like techni- 
color model, based on the Sakai-Sugimoto brane set-up, played an important role in reducing the value 
of the S parameter in that model compared to standard expectations |24j . 
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Note also that we can extract an upper bound, f3u, on the parameter (3 = sin^ a. Namely, 
since the walking region is in the interval p G (1,Pa) and greater values of /3 correspond 
to smaller p\, we can obtain the greates value above which there is no walking, by 



taking pA = 1 in (3.8). This gives: 



/3u = 7.2 X 10" 



(3.9) 



In other words, to have a walking region at all, we need /3 < Clearly, the less stringent 
cutoff arising from |e''^'^ ~ cot^a gives higher upper bound on (3, namely = 3 x 10~^. 



3.2 Computation 



Now we will compute the solutions of (2.27) order by order in small q . We will use a 



different approach compared to [18], which will allow us (due to P <K 1) to calculate 
analytically the constants that remain undetermined by the method used there. Namely, 



integrating (2.27) twice, we obtain the following integral equation: 

rP I 



rp ^ 

C1 + C2 I 77 — rdpi - q 



PO 



Pi 



^' b(fh)^j "*^^2)V^y,A(P2,g^)rfp2jc?Pi, (3.10) 

2 



where ci,2 are arbitrary constants. This can be solved order by order in small q by 
iteration. In particular, in order to find the first order solution, we need to substitute the 
zeroth order tp'^ y in the last term. Note also that, as is easy to realize, under reflection 
z — )• —z of the worldvolume variable z the first term in (3.10) is even, while the second 



one is odd. The third term, on the other hand, is even (odd) when ip^y is even (odd). 
So, to leading order in small g^, we have that: 



1 



pi 



1 



.dpi ( 1 - a{p2j , 

'PO Kpi) V Jpo Jpo KPs) 

where we have taken ci^2 = 1 for easier comparison with [18j as will become clear shortly. 



pi fP2 ^ 

•2) I T^^dp3dp2 ] + 0{q'^ 



(3.11) 



Using that in the walking region the expressions (2.21) acquire the form 

1 



a(p) 

Kp) 



24 ' " V 1 - e'^'-Po-p) ' 



(3.12) 
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where we have introduced the notation cq = ca//?, we can easily compute the integrals in 
(3.11). As a result, we obtain to leading order in small /3: 



S 



and 



2 X 



(3.13) 



(3.14) 



where L = 4(pA — po) + log 4 and so L — 1 > 0. Note that, due to (2.34), we have to 
leading order: Cq ~ -\/3/16. 



Let us make a couple of important remarks about (3.13). First, it clearly implies that 



S* > in our model. And second, because p\ is a function of (5 (recall (3.8)), the overall 
dependence of S on /3 is S' ~ (3^'^/^ / log(/3). As a result, decreasing (3 means increasing S", 
which is rather unexpected. At this point, it is not at all clear whether this is a generic 
prediction of holographic walking technicolor or a strange peculiarity of the model we are 
considering here. We hope to come back to this issue in the future. In any case, we will 
see in the next section that the answer (3.13) is in perfect agreement with the result of 
an independent computation, obtained from summing the contributions of the discrete 
vector and axial-vector modes. 



Finally, for comparison with [18], note that solving the differential equation (2.28) 
iteratively in small one finds: 

^'vAp) = 1 + Ci'^^" + + Cl\') P - ^ (2 - + ^C^2%) + 0{q') , (3.15) 
^^i^J^e~^P. One can extract the coefficients C^2^3 , whose analogues in [18 



a(p). 



where Q 

remained undetermined, by matching (3.15) with (3.11) 

d 
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Then, substituting (3.15) in 



S 



b{p) ^ {dplog^'y - dplog^V) 



(3.17) 



p=pA,q^=o 



we again find (3.13). 
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For completeness, let us record here the answer to leading order in small /?: 



Co 13 



l/2p4po/i _ 



(l-4po) 



48 

^_ co/?i/2e4''o(l-2po) 



CI - 



Co (3 



l/2g4po 



96po 



ct = - 



1 

Po 



ct^ 



12 

Co/3 



l/2g4po 



48po 



(3.16) 
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4 Discrete states 



In this section we will investigate the spectrum of vector and axial-vector mesons. To 
begin with, we will rewrite the relevant equations of motion in a Schrodinger form. The 
change of variables, that enables that, will provide us with a natural worldvolume radial 
coordinate z that runs over the whole D7-D7 embedding. Then we will calculate the 
masses and decay constants of the V and A modes. Finally, we will compute the S 
parameter by summing the discrete state contributions and will find complete agreement 
with the exact result of the previous section, in contrast with the situation in [T6] . 



4.1 Schrodinger form of field equation 



In order to transform the field equation (2.28) into Schrodinger form, let us introduce the 



following variable (very much like the tortoise variable in General Relativity): 



z = ± 



dp , 



(4.1) 



where the "+" sign corresponds to the D7 branch and the "— " sign to the D7 one. 



Substituting (3.12) into (4.1), we find: 



(4.2) 



Note that near p = po the right-hand side of (4.2) behaves as (p — po)^''^- This implies 
that 2^ ~ p — Po for p ~ Po, similarly to the choice of variable z considered in [18]. Let us 
also introduce a function 0(2;) defined via: 



[a(p)6(p)]V^ 



(4.3) 



where [a(p)6(p)]^/^ is viewed as a function of z via (4.2) and, for convenience, we have 



dropped the index n. With these definitions, one can show that (2.26 ) implies the following 
differential equation for <^{z): 



-(^'{z)^V{z)(^{z) =m^(f){z) 



where the Schrodinger potential is given by 



V{z) 



1 z'^-2\ 

'4(^2 + ^)2 



(4.4) 



(4.5) 
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with 

A = ^e'''°/3^/'co . (4.6) 
Since p < Pa, clearly we have that —za < z < za, where 

ZA = T^/J^Ve^M-e^Po ^ J^P'/'e'^^ . (4.7) 



Note also that (3.8) imphes that e^^^ ~ P'^^^ and hence za = 0{p-^/^) » 1. 



In the following subsections we will study the Schrodinger equation (4.4) in order to 
find the mass spectra and decay constants of the vector and axial-vector bosons in our 
model. Before proceeding, let us mention that the above method of transforming equation 



(2.28) to Schrodinger form is equally apphcable to the models studied in [H] and [TB]. In 
their case, one finds that the tortoise variable z has a finite range, i.e. z < z^, for some z^, 
and the Schrodinger potential behaves as V{z) — t- oo for z — )■ z^,. As a result, there is an 
infinite spectrum of vector and axial-vector mesons. 

4.2 Mass spectrum 



Now we will compute the mass spectrum of the solutions of equation (4.4) to leading order 
in small /3. To do that, we will split the range of variation of z (on the positive branch) 
into two subintervals, in which exact solutions can be found, and then we will match those 
two solutions in a region, where they behave in the same way. Namely, we take the two 
subintervals to be (0, z,) and (z,, za), where z, = N\/X with being a number satisfying 
1 « « /3-^/\ 



In the first subinterval, we will keep the form of V{z) in (4.5) unchanged but will 



omit the right-hand side of (4.4). This will be justified a posteriori, since we will see that 



V{z) = 0{(3 ^/^) ^ = 0{(3^^'^) once we compute the mass spectrum. With these 



assumptions, the solution (pi{z) of (4.4) in this interval is: 



(z'^ + A)^/"^, for symmetric (vector) modes, 

0i(z) ~ <( . X (4.8) 

(z^ + A)^/''log ( ' antisymmetric (axial) modes. 

In the second interval, z, < z < za , have that z^ ^ X and so we can approximate 
the potential by V{z) ^ —1/(42;^). Then the general solution 02(2;) of (4.4), satisfying 
the boundary condition 02 (-^a) = 0, is: 

4>2{z)-y/^[Mmz)~RYo{mz)], (4.9) 



16 



where R = Jo{m za) /Yo^m z\) . 

Now, let us match the logarithmic derivatives of (pilz) and 4'2{z) at z = z,, in order 
to obtain the discrete spectrum of vector and axial- vector modes. We start by computing 



from (4.8), to leading order in small \/z,: 



2z, 




\og{2z./^/X) ' 



for vector bosons , 

for axial- vector bosons . 



(4.10) 



Similarly, from (4.9) we find to leading order in small /3: 



2z,- 



2/i 



Mfi) - RYoifi) 



1 + 



2R 



i?[log(m/2) + \og{z,) + 7] - TT / 2 



(4.11) 



where ^ = mz, and 7 is Euler's constant; on the right-hand side of (4.11) we substituted 



the leading small- argument asymptotics of the Bessel functions since /i « 1. Equating 

(4.12) 



(4.10) and (4.11), we obtain the following quantization conditions: 

> Jo{mzA) = 0, 



R = 



and 



R 



TT 

2 



7 + log 



my A 



-1 



for vector bosons 



for axial— vector bosons. 



(4.13) 



Thus, the mass spectrum of vector bosons is given by m^z^ = r;^ with r;^ being the n^^ 
root of Jo- In other words, restoring dimensions, we have: 

2V3 



m. 



V 



Co 



KK 



(4.14) 



where Mkk is a mass scale that is determined by setting = (246 GeV)^ in (3.14). 
Similarly, the mass spectrum of the axial-vector bosons is: 

2v^ 



KK 



(4.15) 



where /i„ is the n root of 



TT 



7 + log 



42a 



(4.16) 



Note that since e'^^^ ~ /3 all of the masses m^'"^ are of order 



Although (4.16) cannot be solved analytically, it is easy to evaluate its solutions nu- 



merically at arbitrary n and zj<^/ y/X = e^*-''^ Through the latter quantity, /i„ is a 
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Figure 1: /i„ (solid lines) and r;{ (dashed lines) for the three lightest vector and axial- vector 
mesons, as a function of e^^P'^~P°^ = / y/X . 



function of the cutoff pA. We argued in Sec. [3] that the cutoff in p should be determined 



by (3.8), instead of the less precise estimate (2.5) given in [IQJ. Despite that, for more 
generality, let us now address how the solutions /i„ depend on the value of p\. The 
choice e^''^/3 /2 = 0.02 is the same as in (3.8), whereas e^P'^P /2 = 1 would correspond 
to Pa = p*. So in order to illustrate the dependence on the choice of cutoff, as well as 
how axial-vector boson masses compare to vector boson masses, we plot p„ and r;[ for 
n = 1,2,3 on Figjl] as functions of e^^''*"''"^ Note that at low orders one has R <K 1. 
Thus the axial wave function is dominated by the term y/z jQ{m z) and the axial- vector 



masses (4.15) are close to the vector ones (4.14). Furthermore, the difference — 
decreases as (5 gets smaller. Nevertheless, the vector bosons remain always lighter than 
the axial ones (at every order), just like in QCD. 



4.3 Decay constants and S-parameter 

In this subsection, we will compute the decay constants gy^ , Qa^ ^he vector and axial- 
vector modes. This will allow us to obtain another estimate for the S-parameter, by 
utilizing the formula [2]: 



n=l 



m 



Unlike [16], we will find that the result from (4.17) agrees completely with the answer 
obtained from the exact method of Section O 
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In HHI, it was shown that 



= / _'^vip)9p[bip)dpipvjdp , 

'D7+D7 

9A^ = I _^A{p)dpHp)dpi)AMp- 

^D7+D7 



(4.18) 



Partially integrating the above, we have: 



Qn 



PO 



PO 



(4.19) 



Now, since dpiljy\q2=o = 0, the vector decay constant gv„ is a total derivative. However, 
dpipWq2=o = ^ 0. Therefore, we find: 

2k 



9a„ 



-2Kb{p) dpilj_ 



pyA„ \p=pA 



+ 



Pa - po 



PA 



Hp) 9pip_ 



P ^Ar^ , 



(4.20) 



PO 



where we have used that ip'^lp^p^ = 1 and 4'a\p=po — 0- Now, partially integrating the 



last term in (4.20) and using that |p=po,pA = 0, we obtain: 

9a„ 



2k b{p) dpijAjp=pA / (dpb) ipA^ . 

Pa - Po Jpo 



PA 



(4.21) 



Since, to leading order in small /3, the walking region has b ^ const, the second term 
vanishes; clearly, though, there will be subleading order corrections arising from that 



term. Thus, the exact expressions (4.18) are equivalent to 



9v,, = -2k6(p)9pV'vJp=pa 
9a„ = -2Kb{p)dpijAjp=pA, 



(4.22) 



that were used in [TB], only in leading order in small /3. 

To compute Qn we need to know the properly normalized wave functions ipn, i-e. 



satisfying (2.25). To find them, we begin with the wavefunctions (pn from Subsection 4.2 



which should be normalized to Sn,m- Note that, in order to normalize 0„ to leading order 
in small /3, we need only consider the region z > z, since the contribution of the region 
z < z, is suppressed by a power of f3. Then we find from (4.9): 



(f^A^z) 



zJo{r^z/zx) 



ZAJi{ri) 
'z{jQ{m^ z) - RYo{m^ z)] 

ZA[Jl{Pn) - RYi{pn)hn 



(4.23) 
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where as before R = Jo{fin) /^o(Atn) and 



7n 



7r/i„[Ji(/i„) - 



(4.24) 



Now one can show that the wavefunctions, normahzed as in (2.25), are given by: 



(4.25) 



Substituting (4.25) into (4.22) we obtain: 



9Vn 



9 Afi 



4 X 3^/^Ke-^P^{ 
4 X 35/^/31/8;. e-^^A(/i„)2 _2 



(4.26) 



Finally, from (4.17) we find for the S parameter: 



S 



K Cn 



33/4 



E 

n=l >- 



7 



-2 



(/^n)2j 



(4.27) 



The sum in this expression can be evaluated numerically to any order n. 

Clearly, as in our earlier discussion regarding fin, the numerical answer for S also 
depends on the choice of cutoff. We found, however, that for every choice of p\ the series 
in (4.27) converges rapidly. The n = 1 term approximates the infinite sum within 20%, 
while the sum of the first four terms, i.e. up to and including n = 4, within 1%. To 
illustrate the cutoff dependence, we plot in Fig|2] the S parameter as a function of (3 
(with Co = v^/16 substituted) for three different values of e'^^''^~''*^ between 0.2 and 0.02. 
Note that there is a perfect agreement between the results in (3.13) and (4.27); the two 
curves coincide within 0.1% for all /3 and px, when we sum over the first 20 states. 

It is worth expressing the various physical quantities directly in terms of F^r. From 



(3.14), we can write: 

Mkk = F^P'/'\^{3co)-'/\ (4.28) 

V K 

which allows us to eliminate Mkk from the expressions for the masses. Namely, substi- 



tuting ( |4.28[ ) into ( |4.14[ ), we obtain: 

V _ „2(p.-pa) 



m 



6L 



K 



16 -5/8q-1/8^J 



C, 







r'^F 



(4.29) 



and similarly for with r:^ ^ pn- Setting = 246 GeV, as appropriate for a technicolor 
model, we have plotted the dependence of the lightest vector and axial-vector masses on 
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0.001 



Figure 2: The S-parameter as a function of (3 at e^'-^* = 0.2 (red), 0.065 (green), and 



0.02 (blue) from (3.13) or (4.27). The black horizontal line represents the experimental 
upper hmit. 



(3 and p\ — p* in Fig|3| Note that so far we have tacitly assumed that Ntf = 1, as in [16]. 
We should point out that increasing the number of flavors would decrease the masses of 



the mesons by a factor of ^/NTp. 

Finally, let us also mention that one can relate directly the physical observables S*, 



and mj^'"^ to each other, without any parameters involved. Indeed, from (|3.14|), (|3.13|) 



(4.14) and (4.15), we obtain 



ra. 



V 




27r 
2^ 




S 



F r 



F-K pn ' 



(4.30) 



where we have used the approximation ^yl — 1/L — )• 1 with L = \og{2z\ / A) ~ — log(/3). 



Note that the expression for is approximate since, due to (4.16), /i„ depends on L as 
well. Setting F^^ = 246 GeV, we obtain for the mass of the lightest technirho boson: 

1 



V 



S 



1.5 TeV. 



(4.31) 



Hence the experimental upper bound for the S-parameter, i.e. S < 0.09 , implies the 
lower bound mY > 5 TeV, for the case of a single flavor. If Ntf > 1, then the lower 
bound on the mass of the lightest vector meson is smaller. We will investigate in more 
detail models with nonabelian flavor symmetry in a future publication. 
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Figure 3: The mass of the hghtest vector (dashed hnes) and axial vector (sohd hnes) 
bosons as a function of (3 at e''*^''^"^*^ = 0.2 (red), 0.065 (green) and 0.02 (blue). 



5 Discussion 

We studied a model of walking technicolor via its dual gravitational description. The 
latter was obtained by embedding probe D7-branes in a type IIB background. We were 
able to compute the electroweak S parameter in two different ways. The first way is via an 
exact expression in terms of non-normalizable modes and the second one is via summing 
the contributions of the discrete vector and axial-vector states. The two answers are in 
perfect agreement, in contrast to previous holographic studies |16] of technicolor models. 
Hence, the issue that was raised there as an explanation for the discrepancy, namely the 
non-decoupling of KK modes, does not occur for us. 

We found an unexpected dependence of the 5* parameter on the length of the walk- 
ing region. It would be very interesting to understand whether this is some kind of a 
characteristic prediction of stringy holographic (as opposed to the AdS/CFT-inspired 
phenomenological) walking models or a peculiarity of the particular model we studied 
here. To address this issue, one would have to try to enlarge the set of examples of 
gravitational duals of walking technicolor. One way to achieve that is to search for em- 
beddings of the techniflavor probe branes, that are different compared to the one of [TH]. 
Such embeddings would have to be more technically involved, but could potentially lead 
to interesting results. Another way to broaden the set of examples is to look for different 
walking backgrounds in string theory. At present, the only known ones are the back- 
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ground we considered here and a slight variant of it [25] (see also [26]). They were found 
at the purely technical level, but at this point there is no good conceptual understanding 
of why they lead to walking. In order to obtain a wider class of solutions, it would be 
rather helpful to achieve a better understanding of those backgrounds. We hope to report 
on that in the future. 

In this work, as well as in [161 H] , the implicit assumption was that there is only one 
pair of probe techniflavor branes. If N^f > 1? then our present results for 5* and would 
be multiplied by a factor of Ntf- On the field theory side, this is just the usual statement 
that each technifermion makes an additive contribution to the spectral functions. On the 
gravity side, that we are studying here, this is due to the fact that for Ntf probe D7-D7 
pairs the total DBI action will be N^f times the action in (2.19). The presence of this 



factor on the right-hand side of (3.14) would imply a decrease of the discrete spectrum 
masses by \/Ntf- More importantly, in addition to the inclusion of an overall multiplier, 
the proper consideration of the Ntf > 1 case requires taking into account the nonabelian 
nature of the techniflavor field strength. This is crucial in order to be able to compute 
the T and U electroweak observables, which we will address in a future publication. 

The result we found for the S-parameter is positive-definite. However, although we 
can decrease it by varying the parameters of our model, we cannot reach S = within our 
approximations. Despite that, our result is below the lower bound that was conjectured 
in [27] based on the existence of an IR fixed point, the Banks-Zaks point [2H]- One 
should keep in mind, though, that our model is, quite likely, in a different universality 
class compared to the more familiar IR-fixed-point ones [3]. A clear indication for this 
is provided by the fundamentally different role the number of technifiavours is playing in 
the two cases. Namely, in the IR fixed point considerations, it is crucial for the number of 
technifiavours to be comparable to the number of technicolors (i.e. Ntf ~ Ntc), in order 
to achieve a nearly conformal regime. In our considerations, on the other hand, one always 
has Ntf ^ Ntc- Hence, our gravitational description is not a dual of a standard WTC 
model [3], but instead a dual of a novel kind of walking model. As a consequence, our 
positive-definite answer cannot be viewed as evidence against the possibility, considered in 
[29118]. of having models with negative S parameter. In fact, it would be rather interesting 
to explore whether one could obtain 5* < from a different kind of gravitational dual, 
probably with A^^^ ~ Ntc- 

In order to reproduce on the gravity side the relation Ntf ~ Ntc, used in the IR 
fixed point models, one would have to go beyond the techniflavour probe approximation. 
The standard trick to construct a background with backreacted technfifiavour branes is 
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to smear the latter along the transverse directions, in order to simplify the supergravity 
field equations. However, such a smearing is incompatible with the U-shaped embedding 
that realizes geometrically chiral symmetry breaking, which was of crucial importance for 
us. Hence, at present, finding a backreacted gravitational dual of walking technicolor, in 
which Ntf is of the order of Ntc, is a technically-challenging open problem. 
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